FERMILAB-CONF-04-196-E 
February 7, 2008 



Low-Energy Neutrino Beams with an 
Upgraded Fermilab Proton Driver 



S.J. Brice, S. Geer 

Fermilab 

K. Paul 

University of Illinois 

R. Tayloe 

Indiana University 



Abstract 

The beam properties of low energy neutrino beams that would be possible at an 
upgraded 2 MW Fermilab 8 GeV Proton Source are described. In particular three 
options are considered; (i) a MiniBooNE-like beam using a conventional neutrino horn, 
(ii) a decay at rest neutrino source, and (iii) a neutrino beam downstream of a solenoid 
decay channel. In all three cases the fluxes are sufficiently large to provide an interesting 
physics program. Some physics examples are considered. 
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Figure 1: The simulated and v e flux energy spectra at the MiniBooNE detector. (From 
0)- 

1 Introduction 

The Fermilab Long Range Planning Committee P has considered which facilities and as- 
sociated physics programs might provide an exciting future for the laboratory beyond the 
Tevatron Collider era. Two accelerator-facility-based physics programs have emerged as 
leading candidates for Fermilab's future: (i) a Linear Collider and (ii) an expanded neutrino 
oscillation physics program driven by a new high intensity 8 GeV proton source delivering a 

2 MW beam on target [2j. A small fraction of the beam from the new proton source would 
be injected into the Fermilab MI to provide a higher energy (up to 120 GeV) proton beam, 
also with 2 MW beam power on target. This scheme would enable the intensity of the NuMI 
neutrino beam to be increased by a factor of about 5 for on- axis and/or off- axis long-baseline 
experiments, whilst at the same time providing an almost 2 MW proton beam at 8 GeV for 
lower energy neutrino experiments. In this note we describe the properties of various low 
energy neutrino beams that could be created with a 2 MW proton beam at 8 GeV. 

2 A High Intensity MiniBooNE-Like Beam 

The existing MiniBooNE beamline and detector [Hj provide a useful point of comparison 
for neutrino beams at a future proton driver. MiniBooNE takes 8 GeV protons from the 
FNAL Booster, steers them into a Beryllium target, and uses a magnetic horn to focus the 
resulting mesons. The mesons are given 50m to decay before hitting a concrete absorber. 
The resulting neutrino beam crosses 500m of dirt before reaching the MiniBooNE detector, 
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Table 1: Yearly rate of event types in a MiniBooNE-like detector at an 8 GeV, 2 MW 
proton driver. (Data from 

a 6.1m radius spherical tank of mineral oil viewed by 1520 PMTs, looking at the Cerenkov 
and scintillation light created when neutrinos interact and produce charged particles in the 
tank. Fig. ^ shows the simulated and v e energy spectra at the MiniBooNE detector. 

MiniBooNE's design proton intensity calls for 5 x 10 12 protons per spill, repeating at 5 
Hz. For a 10 7 -s operational year this corresponds to 2.5 x 10 20 P/yr. An 8 GeV, 2 MW 
proton driver would deliver 1.6 x 10 22 P/yr. Therefore, the event rates with a 2 MW proton 
driver can be obtained by scaling up the MiniBooNE event rates by a factor of 160/2.5 = 64. 
The rates for the dominant neutrino interaction channels are shown in Table H 

3 Decay-at-Rest Neutrino Beams 

The principal advantage of a neutrino beam from muon and pion decay at rest over beams 
produced in other ways is the extremely well known neutrino energy spectrum. The main 
disadvantage is that the flux is isotropic and drops as 1 /R 2 compared to the high degree of 
collimation of other beam types. 

3.1 The Beam 

A multi-GeV proton beam will interact in a target to produce pions and kaons. If the target 
is sufficiently large the charged pions will lose energy and come to rest. As they slow down 
some of the pions will decay to produce daughter muons, which will also lose energy and 
come to rest. The fate of the charged pions and muons will depend upon whether they are 
positively or negatively charged. The stopped positively charged pions and muons will decay 
at rest, producing a "prompt" flux from 7r + decay followed by a "delayed" v e and 
flux from the much slower u + decays. This time structure will be evident if the incoming 
primary proton beam bunch is short compared to the muon lifetime. In contrast, before 
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Figure 2: The normalized energy spectra of the neutrinos produced from ll + decay at rest, 
where x = is the fraction of the maximum available energy given to the neutrino during 
decay. 



they can decay most of the stopped negatively charged pions and muons will be absorbed 
by the atomic nuclei (the relevant capture rates are much larger than the decay rates). The 
suppression of negatively charged pion decays is enhanced in high-Z target materials for 
which the capture time is very short. 

Pion and muon decays at rest (DAR) produce an isotropic flux of low energy u e , u^, 
and i>fj, with energies between a few MeV and half the meson rest-mass energy. Neutrino 
spectra from DAR are shown in Fig. |2l Since the flux is isotropic, the challenge is to obtain 
sufficiently high intensities to be of interest to cutting-edge experiments. The key to meeting 
this challenge is to start with a very high intensity proton source. For this reason there has 
been some interest [3] in the possibility of exploiting the soon-to-be-completed Oak Ridge 
SNS facility as a low energy neutrino source. The SNS will provide a 1.4 MW primary proton 
beam at 1.3 GeV with a pulse length of 695 ns. The pulse length is short compared to the 
2 /as muon lifetime, and hence suitable time cuts can separate the neutrino flux produced 
from the decaying pions from the flux produced by decaying muons. During the duration of 
the proton pulse (from t = ns to t = 695 ns) 96% of the positive pions decay. From t = 695 
ns to t = 5 lis the remaining 4% of the positive pions decay and 74% of the muons decay. 

To a first approximation, at fixed beam power, low energy pion production rates are 
independent of proton beam energy. Hence a first guess at the neutrino fluxes from a DAR 
source at a 2 MW Fermilab proton driver can be based on DAR studies for the SNS, scaling 
upwards by a factor of about 1.4. To obtain a better estimate we have used the MARS code 
and simulated particle production, interaction, energy loss, and decay for an 8 GeV primary 
proton beam incident on a carbon target with dimensions sufficiently large to contain all of 
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Table 2: Neutrino production rates from a carbon target and the event rates in a MiniBooNE- 
like detector placed 90 degrees off axis, 60 m from the target. The production rates are 
computed using MARS with an 8 GeV, 2 MW proton driver for a 10 7 -s operational year. 
The event rates in the detector are found by scaling from G. VanDalen's paper j3] with 
fluxes computed by MARS for an 8 GeV, 2 MW proton driver in operation for a 10 7 -s year. 
Dependence of the rates on the neutrino reconstruction efficiency, e rec , is explicitly indicated. 
The symbol v indicates all neutrino and antineutrino flavors. 



the produced pions. For comparison we have also simulated a 5 GeV proton beam incident 
on the target. To calculate neutrino event rates we have used the calculated SNS DAR rates 
and scaled these numbers to correspond to our calculated fluxes, which assume a 2 MW 
beam at 8 GeV (or 1.25 MW at 5 GeV), and a MiniBooNE-like detector placed 90 degrees 
off proton-beam axis at a distance of 60 m from the target. The resulting event rates are 
shown in Table |2]for an operational year (10 7 s). Note that the 2 MW, 8 GeV proton source 
produces DAR neutrino fluxes that are a factor of 2.5 higher than a potential SNS DAR 
source, somewhat better than suggested by our naive estimate. 

3.2 The Physics 

There are a number of fundamental physics measurements one can consider with a decay at 
rest beam. If MiniBooNE should confirm the indication of neutrino oscillations reported by 
the LSND experiment then a new generation of experiments investigating the oscilllation in 
fine detail and, in particular, its energy shape will be needed. Besides having an extremely 
well characterised neutrino spectrum the backgrounds in an LNSD- or MiniBooNE-style 
detector would be extremely low from a proton driver DAR beam (the beam related back- 
grounds are extremely low and a small machine duty factor would ensure very low cosmic-ray 
related backgrounds). This arrangement might be the ideal followup to MiniBooNE if that 
experiment should see an oscillation signal. Using the rates from Table El a MiniBooNE sized 
detector operating at an 8 GeV proton driver would see many hundreds of oscillations events 
per year (varying with the exact oscillation parameters) with about 5 beam unrelated back- 
ground events and a couple of beam related backgrounds (including muon misIDs). These 
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numbers are scaled from The low backgrounds and well-known spectrum at a DAR source 
more than make up for the drop in neutrino cross-section at the resulting comparatively low 
neutrino energies. 

It has been suggested [2j that a neutrino beam from muon DAR can be used to search for 
a non-zero #13, the only unmeasured angle in the three neutrino mixing matrix. The effect 
of 0i3 on v e appearance is maximised when the oscillation baseline is at the atmospheric 
Am 2 . For muon DAR energies this puts the detector about 10km from the source and a 
very massive detector (on the scale of the proposed HyperK) would be required. Whilst 
this may be possible it is unlikely to happen on a timescale competitive with reactor or 
superbeam experiments. 

Table|2]shows a healthy rate of neutrino electron scattering in a MiniBooNE scale detector 
at an upgraded FNAL proton driver. These events can be used to make a very low Q 2 
measurement of sin 2 6^, the Weak mixing angle, by measuring the cross-section ratio 

R = , "l V - e) , , . (1) 

Neutrino electron scattering events can be identified by looking for very forward going elec- 
trons in the detector, and events from the numerator of the R ratio come in time with the 
beam whereas events in the denominator come later (in general). Assuming all neutrinos 
have the same flux shape, and neglecting terms of order m e /E u , this ratio calculated in the 
SM is 

0.75 -3 sin 2 fl w + 4sin 4 6 W 
1 + 2 sin 2 6 W + 8 sin 4 9 W ' ^ ' 

Based on the work of |B] a measurement of R in this way at an upgraded FNAL proton 
driver would yield a 1-2% measurement of sin 2 ^. This measurement would be at low Q 2 
compared to other weak mixing angle measurements. In fact the only other measurement in 
this Q 2 range comes from the totally different regime of atomic parity violation. In view of 
the recent indications of non-SM behaviour in the extraction of sin 2 6 W from neutrino electron 
scattering at higher Q 2 jH] this measurement from a decay at rest beam from an upgraded 
FNAL proton driver is very attractive. 

Neutrino electron scattering can also be used to search for a neutrino magnetic moment 
as it adds an electromagnetic component to the otherwise weak neutrino electron scattering 
cross-section. This component gets stronger with decreasing electron recoil energy and so 
low energy neutrino sources such as reactors or stopped muon decay are the best places to 
go looking for it. The well known energy spectrum of decay at rest neutrinos is a particular 
advantage when looking for a distortion in the recoil electron energy spectrum caused by a 
neutrino magnetic moment. 

4 Neutrino Beam from a Solenoid Decay Channel 
4.1 The Beam 

Neutrino beams are usually produced using a magnetic horn to sign select and focus the 
secondary charged pions into an approximately parallel beam that then propagates within 
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Figure 3: This figure depicts the typical layout of a solenoid-based pion capture and decay 
channel used for generating neutrino beams. 

a long large aperture decay pipe. An alternative approach is to use the optics that has 
been proposed for the upstream stage of a neutrino factory, in which a pion production 
target is within a high-field solenoid with a field strength Bq, sufficiently high to radially 
contain essentially all of the produced charged pions. This means that, if px is the maximum 
transverse momentum of the produced pions, then B must be chosen such that 



where Rq is the radius of the beampipe around the target. Note that for an 8 GeV primary 
proton beam incident on a carbon target the secondary pion transverse momentum peaks at 
about 150 MeV/c. For the design we have explored Rq = 7.5 cm and Bq = 20 T. This system 
will capture particles with px < 225 MeV/c, resulting in a captured beam with significantly 
broader energy spectrum than one captured with a horn. 

It would be unnecessarily expensive to extend the 20 T solenoid for the full length of the 
decay channel downstream of the target. Therefore, downstream of the target the charged 
pions, and their daughter muons, are confined radially within a lower field (B = 1.25 T) 
larger radius (R = 30 cm) solenoid decay channel. To ensure efficient transfer of the beam 
from the high-field to the low-field solenoids, the distance over which the field is decreased 
must be approximately 10 m or more. Within this section the focusing field strength B(s) is 
adiabatically decreased whilst the beampipe radius R(s) is increased, keeping the magnetic 
flux within the beampipe, <J> oc B(s)R(s) 2 , constant. The continual transverse focusing of 
the beam throughout the entire length of the decay channel allows for the marked increase 
in muon yields over those from horn designs. A typical layout is shown in Figure E3 

The charged pions within the solenoid channel decay via 7r + — > yU + v ^ and 7r~ — > yT u^, 
with a lifetime of r n = 26 ns. The daughter muons decay with a lifetime of r M = 2.2 fis via 
/i + — > e + u e and fi~ — > e~ v e v^. Due to the significant difference between the pion and 
muon lifetimes, the muon content of the beam will steadily increase over the first ~75 m 
downstream from the target. At greater distances the muon content of the beam decreases 
as the muons decay. The muon and charged pion content of the beam is shown in Table El as 
a function of distance downstream from the target, s, for an 8 GeV, 2 MW primary proton 
beam incident on a carbon target with a 25 m adiabatic section following the 20 T capture 
solenoid (Figure EJ). The corresponding evolution of the charged pion and muon spectra is 
shown in Figure 0] Note that the simulated layout does not sign select the pion and / or muon 
beam. Hence, downstream of the decay channel, the neutrino beam contains neutrinos from 
7r + ,7r~,/i + and \i~ decays. However, the design can be modified by placing suitable large 
angle bends at any desired positions along the decay channel. Thus sign selection bends, 
and bends to enhance of the muon DAR component, are possible options. 
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Charged Particle Fluxes in the Decay Channel with 1.6 x 10 22 POT at 8 GeV 
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Figure 4: Charged particle fluxes at various positions down the 100 m decay channel with a 
30 cm beampipe radius. We assume an 8 GeV, 2 MW proton driver in operation for a 10 7 -s 
year. 



7 



Neutrino Fluxes in the Decay Channel with 1.6 x 10 POT at 8 GeV 
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Figure 5: Neutrino fluxes at their point of origin in the four 25m segments down the 100 m 
decay channel with a 30 cm beampipe radius. We assume an 8 GeV, 2 MW proton driver 
in operation for a 10 7 -s year. 
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The flavor content of the neutrino beam depends upon the length of the decay channel 
(Tableland Figure|3J). Early in the channel pion decays dominate. In the last 25 m segment 
of the decay channel the pions constitute a few percent of the charged particle beam, and the 
number of muon and pion decays are similar. Figures [3 and El show the neutrino fluxes from 
neutrinos born in the last 50 m of the decay channel at a distance of 100 m from the end of 
the decay channel. The fluxes are averaged over a 10 m radius from the axis of the beam. 
Figure [7| shows the neutrino fluxes produced solely from pion decays, and Figure |H1 shows the 
corresponding neutrino fluxes from muon decays. It is worth noting that, while the number 
of pion and muon decays in the last 50 m of decay channel are comparable, the neutrinos 
coming from pion decays are considerably more collimated because of their comparatively 
higher energy in the last 50 m of the decay channel (Figures El and |UJ). 

To calculate event rates we consider a MiniBooNE-like detector placed on axis at a 
distance of 100 m from the end of the decay channel. The rates for different interactions and 
neutrino types, assuming an 8 GeV, 2 MW proton driver in operation for a 10 7 -s year, are 
shown in Table El The cross-sections on CH 2 are computed using the NUANCE code jU] . 

4.2 The Physics 

As can be seen in Figure the iy^) flux from tt + (ir~) decays-in-flight resulting from this 
neutrino beam are large. The resulting neutrino event rates in a MiniBooNE-like detector 
100m from the end of the decay channel will be unprecedented and could be used to precisely 
measure the neutrino oscillation parameters for both v and v in the L/E ~ 1 range as is 
currently being investigated by the MiniBooNE experiment [3]. 

If MiniBooNE confirms the LSND result, this will indicate neutrino oscillations in the 
L/E m 1, Am 2 m leV 2 range, and a systematic study of this effect, for both v and u, will be 
required. A solenoid-based neutrino beam with a large angle bend to sign select the pions 
would be an excellent source for these studies. The bend will allow for a relatively pure v or 
v beam, depending on the polarity of the bending magnet. The polarity reversal could take 
place on the timescale of hours or minutes, thus reducing systematic uncertainties between 
v and v running arising from detector response. 
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Table 3: The number of charged particles per incident proton in the beam at various distances 
downstream from the target. These numbers are computed by MARS for an 8 GeV, 2 MW 
proton driver such as the proposed Fermilab upgrade. 



9 



Neutrino Radial Distribution at the Detector with 1.6 x 10 22 POT at 8 GeV 
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Figure 6: This figure shows the fluxes, at a distance of 100 m from the end of the solenoid 
decay channel and averaged over lm-radius annuli, of the neutrinos produced in the last 50 
m of decay channel. We assume an 8 GeV, 2 MW proton driver with operation over a 10 7 -s 
year. 
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Neutrino Fluxes from 7r ± Decays at the Detector with 1.6 x 10 POT at 8 GeV 
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Figure 7: This figure shows the fluxes, at a distance of 100 m from the end of the solenoid 
decay channel and averaged over a radius of 10 m, of the neutrinos produced solely from 
pion decay in the last 50 m of decay channel. We assume an 8 GeV, 2 MW proton driver 
with operation over a 10 7 -s year. 
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Neutrino Fluxes from ^ Decays at the Detector with 1.6 x 10 22 POT at 8 GeV 
T 




2 

Energy (GeV) 

Figure 8: This figure shows the fluxes, at a distance of 100 m from the end of the solenoid 
decay channel and averaged over a radius of 10 m, of the neutrinos produced solely from 
muon decay in the last 50 m of decay channel. We assume an 8 GeV, 2 MW proton driver 
with operation over a 10 7 -s year. 
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Table 4: The number of neutrinos produced per incident proton in the four 25 m segments of 
the 100 m decay channel, as well as those neutrinos produced that reach the detector. These 
numbers are computed by MARS for an 8 GeV, 2 MW proton driver such as the proposed 
Fermilab upgrade. 



For these experiments, two large (~ lkton) detectors, would be required at distances on 
the order of 100m and 250m (the exact locations would be determined by the neutrino data 
at the time of experiment construction). In addition, a relatively small detector (« 20 tons) 
would be situated at a near location (« 25 — 50m). This detector would be highly segmented 
to precisely monitor the flux x cross section and would have charge identification to measure 
the "wrong-sign" component of the neutrino beam. 

With this configuration of neutrino beam and detectors, simultaneous precision measure- 
ments of Up — > v e , — > i> e , — > v x , and — > v x in the Am 2 =.1-1 eV 2 region would 
be achievable. With the correct arrangement of detectors, an oscillation measurement using 
neutrinos from /i + that decay-at-rest in the beam stop could also be conducted in parallel. 
These neutrinos have a different time structure due to the long-lifetime of the /i + and are of 
much lower energy as described above. These factors will allow for a clean separation of the 
neutrino sources. 

In addition to servicing the oscillation measurements, the near detector would be a rich 
source of non-oscillation neutrino physics. A 20 ton detector at 25m would collect events 
at 10% of the rate shown in Table 5. This would allow measurements of neutral- and 
charged- current processes with unprecedented precision. These details of these processes 
(cross-sections, energy- and Q 2 -dependence) are needed for oscillation measurements and 
provide valuable information on nucleon structure. 

5 Conclusions 

We have considered three options for generating a low energy neutrino beam using an up- 
graded 2 MW Fermilab Proton Source; (i) a MiniBooNE-like beam using a conventional 
neutrino horn, (ii) a decay at rest neutrino source, and (iii) a neutrino beam downstream 
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Table 5: The above table shows the total integrate events rates from neutrinos produced 
in the last 50 m of the decay channel from pion and muon decays. The MiniBooNE-like 
detector is placed at 100 m from the end of the decay channel. The fluxes are computed by 
MARS for an 8 GeV, 2 MW proton driver in operation for a 10 7 -s year and the cross-sections 
on CH 2 are computed by the NUANCE code 0. 



of a solenoid decay channel. In all three cases the fluxes are sufficiently large to provide an 
interesting physics program. The upgraded MiniBooNE-like beam would be a factor of 64 
more intense than the existing beam. The decay at rest beam would offer an intensity that 
is a factor of 2.5 higher than the corresponding intensity at an equivalent source at the SNS. 
The solenoid decay channel beam has some characteristics which would make it particularly 
interesting if MiniBooNE were to confirm the LSND result, provided sign-selection can be 
implemented efficiently in a (bent) solenoid. 
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